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TECBNICAL NOTE 2571 


APPLICATION OF THE VON KARMAN MOMENTUM THEOREM 
TO TURBULENT BOUNDARY LAYERS 
By Jerold M. Bidwell 


SUMMARY 


A study is made of the Von Kannan momentum theorem with respect to 
its application to turbulent boundary layers in a positive pressure 
gradient. Although the Von Karman momentum theorem for turbulent boundary 
layers contains momentum terms due to the fluctuating motion as well as 
momentum terms due to the mean motion, the general practice has been to 
neglect the momentum terms due to the fluctuating motion. Data were 
obtained from Schuhauer and Klehanoff (NACA TN 2133) and Ludwieg and 
Tillmann (NACA TM I 285 ) with which both the terms due to the mean flow 
and the terms due to the fluctuating flow could he evaluated. The results 
indicate that the streamwise derivative of the turh\LLent longitudinal 
momentum puTuT may he large near separation and therefore should he 
considered when the Von Harman momentum theorem is used for turbulent 
boundary layers near separation. 


INTRODUCTION 


Because the Von Karman momentum theorem (reference 1) is a relation 
between the changes in boundary -layer momentum and the external stresses 
of pressxire gradient and wall shear, it applies equally well for 
turbulent boundary layers as for laminar boundary layers. For turbulent 
boundary layers, however, problems arise in the interpretation of 
momentum and in the application of the theorem to experimental data. 

In a turbulent flow, the net motion is a combination of molecular 
motion, turbulent motion, and stream motion. Reynolds (reference 2) 
showed that the periods of oscillation associated with these three 
phases of motion are distinguishable by means of appropriate time averages. 
In aerodynamics, the primary interest in the molecular motion lies in 
its integrated effects, that is, pressure, temperature, and viscosity. 

In the study of turbulent motion, the fluid is generally considered to 
he a continuum and the net velocity at any point is assumed to he resolved 
into a fluctuating component and a mean component which is independent 
of time. 
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If the velocities appeared in the equations of motion only to the 
first power, a long-time average of the fluctuating component would he 
zero and the equations would contain terms involving only the mean 
velocities. Because, however, the momentum involves the squares of the 
velocities, momentum terms involving time averages of the squares and 
products of fluctuating components must he considered inasmuch as they 
are not zero. The expression for the Von Kannan momentum theorem as 
applied to turbulent boundary layers therefore contains momentim terms 
for the mean flow and momentiam terms involving the fluctuating components. 

Methods (references 3 to 5 ) for calculating the development of a 
turbulent boundary layer have, in general, used an empirical equation 
for wall shearing stress, a known pressure distribution, an en^jiriceil 
equation relating the shape of the velocity distribution to the external 
forces, and the Von Karaan momentum theorem to give the changes in 
boundary-layer momentum. The problem is determinate if the changes in 
boundary-layer momentum are assumed to be ascribed to the mean velocity 
alone. Calculations made with these assumptions have given good results 
in many instances but recently discrepancies have been noted, especially 
in the region of separation. 

Discrepancies have also been noted when the Von Ke^m^n momentum 
theorem was used to obtain wall shearing stress (references 6 to 9) • 

By assuming that the changes in boundary-layer momentum could be 
completely described by the mean flow, measurements were made of the 
mean velocity and pressure gradient, and the Von Kannan momentum theorem 
was used as a balance to give the wall shearing stress. This method 
gave values of wall shearing stress that increased in a positive pressure 
gradient which is in contradiction to the data of references 10 to 12 
which indicate that the wall shearing stress decreases in a positive 
pressure gradient. 


The neglect in the Von Karman momentum theorem, as it is usually 
applied to turbulent boundary layers, of the momentum terms involving 
the averages of the fluctuating velocities is a possible explanation 
of the observed discrepancies. This possibility is examined in the 
present paper. With the use of the Kavier-Stokes equations as a starting 
point, an integral relation that Includes terms resulting from the 
fluctuating motion as well as those associated with the mean flow is 
derived for the turbulent boundary layer. With the use of this relation 
and the experimental data given by Dryden (reference 10) and Schubauer 
and Klebanoff (reference 11), together with measurements of the wall 
shearing stress made Independently by Ludwleg and Tillmann (reference 12) , 
an attempt is made to evaluate the relative order of magnitude of the 
terms resulting from the fluctuating motion as compared with those 
resulting from the mean motion. 
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After the analysis presented herein had been completed, research 
reported by Wallis (reference I3) in Australia was made available. 

Wallis surveyed the equations of energy for turbulent boundary layers 
and suggested that "where the turbulent boundary le^er flow is being 
rapidly accelerated or decelerated, the von Karman momentum equation will 
be inaccurate as it does not take account of the internal process by 
which the energy of mean motion is altered Independently of an external 
force such as skin friction." Reference is made by Wallis to a private 
communication from B. G. Newman in which the problem is approached from 
momentum rather than energy considerations. Newman STiggests that the 
Reynolds normal stresses modify the Von Karman momentum theorem near 
separation. Apparently, data were not available with which a suitable 
check of the hypothesis could be made. 

An analysis similar to that presented herein but much less detailed 
is included in a recent paper by Rubert and Persh (reference l 4 ) . 


SYMBOLS 

H boundary -layer -velocity-profile shape parameter (S*/ 0 ) 

p absolute static pressure at any point 

t time 

u streamwise velocity at a point in boundary layer 

U streamwise velocity at outer limit of boundary layer 

V normal velocity at a point in boundary layer 

X streamwise coordinate 

y normal coordinate 

5 boundary -layer thickness 

5 * displacement thickness 

0 momentum thickness 

H coefficient of viscosity 

V coefficient of kinematic viscosity (ki/p) 
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p fluid density 

T viscous shearing stress 

Subscripts: 

m reference condition 

0 wall boundary 

6 outer limit of boundary layer 

Primes denote fluctuating quantities and bars denote time averages. 


AHAIYSIS 


The Wavier -Stokes equations for a two-dimensional flow in which 
the viscosity density are assimied constant and the body forces are 
negligible are 




+ 



V 


Su 

By 


1 ^ 
P Bx 




|l+u^ + v^=-i^ + ^ 

Bt By p By 


( 1 ) 

( 2 ) 


The associated equation of continuity is 


Bu ^ ^ 
Bx By 


= 0 


( 3 ) 


By resolving the instantaneous values into mean and fluctuating components, 
by considering the mean flow steady and taking a long-time average of 
the fluctuating components, by considering the Prandtl boundary -layer 
assumptions valid for the mean flow up to separation, and by using the 
continuity equation, equations (l) and (2) may be written as 




i. ^ + V 
J5" Bx ^ 


+ 


(4) 
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Su'v' dv*Y' _ 1 ^ 

^ ^ P Sy 

In the immediate neighborhood of separation, the validity of the houndary- 
layer assumptions may he open to question. The extent of the error 
introduced hy making the houndary-layer assumptions under such circum- 
stances is, however, not known. In any case, hy retaining the tenns 
due to turbulence, it is not Implied that these terms are necessarily 
larger than those that have heen neglected hut that the magnitude of 
the effect of the terms due to turbulence is to he Investigated. By 
combining equations (4) and (5) and hy integrating through the boundary 
layer, the following equation of momentum is obtained (see appendix): 



U2 



Sx^ 


«ay 


( 6 ) 


The Von Karman momentum theorem for laminar boundary layers is 


tq ^ 
pu2 ' <3x 


(H + 2) 


e ^ 

U dx 


0 


The usual practice has heen to consider this equation equally valid for 
turbulent flow. 

Whether the Von Karman momentum theorem can he used for turbulent 
boundary layers hy considering only the momentums of the mean flow can 
he ascertained hy experimenteil evalxiation of the momentum terms on the 
right side of equation ( 6 ). 


PRESENTATION OF DATA AND EVALUATION OF THE MOMENTUM TEEMS 


/ / 

Data are presented from which the terms in the Von Karman momentum 
theorem can he evaluated and the order of magnitude of t}ie terms due to 



6 


NACA TN 2571 


turbulence can be compared •with the magnitude of the terms due to the 
mean flow. 


Mean-Flow Data 

Presentation of data and discussion of accuracy . - The mean-flow 
data of Schubauer and Klebanoff (reference 11) have been plotted as 
(U/Um)^, 0, and H as fimctlons of x in figures 1, 2., and 3} 

respectively. The original velocity data were studied in some detail, 
and the scatter of data is believed to be a good indication of the 
random errors in the measurements and in the integration procedinre. 

A systematic error occurs in the boundary-layer data because, in a 
turbulent air stream, pitot-tube measurements give values of velocity 
which are in excess of the mean local velocity. This eiror has been 
discussed by Goldstein in reference 15 and is shown to be 



The velocity data given by Schubauer and Klebanoff have not been 
corrected for the turbulence but sufficient data are presented to permit 
this computation. The magnitude of this correction to the velocity 
is illustrated in figure U for a point just before separation 
(x = 25 .^ ft). The difference between the curves is somewhat greater 
than differences that might be Indicated by random errors alone. 

The velocity profiles have been corrected and integrated and the 
corrected values of 0 are plotted in figure 2. At small values of x, 
that is, less than 21.5 feet, the correction to the velocity profile 
was negligible. Because the correction resulted in a decrease of 
velocity, the displacement thickness was increased and the values of 
corrected H were considerably higher. The corrected values of H 
are plotted in figure 3 aud refaired with a dashed line. The corrected 
veilues of these parameters have been used in the subsequent computations. 

Comparison of wall shearing stress as obtained by different methods . - 
The values of wall shearing stress estimated by Schubauer and Klebanoff 
are compared in figure 5 with the values of wall shearing stress obtained 
by use of the formulas of Squire and Young (reference I 6 ), Falkner 
(reference 17) ^ aud Ludwieg and Tillman (reference 12). The flat-plate 
(zero pressure gradient) formula of Squire and Yoxing has been used 
extensively in the United States in conjunction with the method of 
Von Doenhoff and Tetervln (reference 4) and the flat-plate fonuula of 
Falkner has been used in England in conjunction with the method of 
Garner (reference 5)- The formxila of Ludwieg and Tillman was the result 
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of some experiments made to determine wall shearing stress in a pressure 
gradient with the aid of a calibrated heat-transfer instrument. 

Although the estimated values of wall shearing stress given by 
Schubauer and KLebanoff became zero at separation, at lower values of x 
these values of wall shearing stress are considerably Mgher than those 
given by the flat-plate formulas. A positive pressure gradient tends 
to decrease the velocity in the boundary layer and hence to decrease 
the normaJ. gradient of velocity at the wall. Because the flow in the 
immediate neighborhood of the wall is laminar, the decrease in normal 
velocity gradient at the wall is directly proportional to the decrease 
in wall shearing stress. For ail values of x larger than 17* 5^ 
pressure gradient is positive; so the wall shearing stress would be 
expected to be less than that obtained by flat-plate formulas. Schubauer 
and Klebanoff (reference 11) estimated these values of wal l shearing 
stress by fairing the pu'v' data to the region of the wall and then 
extrapolating the wall shearing stress by continuing the fairing to the 
wall with a slope equal to the pressure gradient. Although the trend of 
the data obtained by this method is correct, inaccuracies in the extra- 
polation resulted in absolute values that are apparently too large. 

The distribution of wall shearing stress obtained by using the 
fonnula of Ludwieg and TlUmann appears to be the most reasonable from 
consideration of the experimental method used and from comparisons with 
well-established results for the zero-pressure-gradient case. Inspection 
of the Ludwieg and Tillmann formula 


tq _ 0.123 1 

py 2 - ^q0.678h j^O.268 

(where Rg is the Reynolds number based on O) reveals that H must 
be infinite for the wall shearing stress to disappear and an infinite H 
is, of course, physically unreasonable. At separation, however, the 
indicated value of wall shearing stress is small enough to be considered 
negligible . 

Substitution of mean-flow data in the Von Karman momentum theorem . - 
Values of -d0/dx and (H + 2) ^ ^ were obtained by use of figures 1, 

2, and 3 and are plotted as a function of x in figure 6. The wall 

shearing stress tqJ pU^ computed from the formula of Ludwieg and Tillmann 
is also plotted in figure 6. If in the application of the Von Karmen 
momentum theorem to turbulent boundary layers only mean-flow quantities 
need be considered, the sum of the terms would be zero. The sum of the 
mean-flow terms is plotted as a function of x in figure J, where it 
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is shown that for a large range of x consideration of mean-flow 
quantities alone gives good, restilts; near separation, however, the sum 
has a large valxie. Good, results are obtained, up to a value of H of 
approximately 1-5 • 

Because changes in the fairing of 6 will alter the value of 

it is of interest to know how much different the values of 0 woiild 

have to he to make the sum of the mean-flow terms zero. The — required 

dx 

to make the s\mi of the mean-flow terms zero was integrated, with respect 
to X from 17 • 5 feet to 25.^ feet and. the results are plotted, as a 
d.ashed. line in figure 2. The differences between the experimental curve 
and the required curve could not be attributed to differences in fairing. 


Pluctuating-Flow Data 


Presentation of data.- Values of ^ reported by Schubauer and 

^ 

Klebanoff have been plotted as a function of x and y in figure 8. 

The data were plotted "three dimensionally" to facilitate fairing. 
Successive planes of x = Constant were indicated along a axis at 
distances equal to particular values of x. By this method of plotting, 
curves of y = Constant can be used to check the fairing and to aid 
the visualization of the surface. This method of plotting and fairing 
can be xised because the function indicated by the fluctuating- flow data 
is regular up to separation. 

Evaluation of the momentum terms due to turbulence .- The faired 

data of figure 8 and figure 1 were used to obtain the values necessary 
to evaluate the term 


u2 Jo 

which is plotted as a function of x in figure 7- The magnitude of 
this term is affected only moderately by the fairing of the data. The 
value of this function is small at lower values of x and becomes large 
only in the region of separation. A comparison of the value of this 
function with the sum of the mean-flow terms obtained from figure 6 

Indicates that, at least qualitatively, the term — / dy, together 

TJ2 Jq Sx 

with the mean-flow terms may be sufficient to describe the momentum 
balance in the turbulent boundary layer up to separation. A more precise 



2P 


MCA TW 2571 


9 


discussion of the in^jortance and absolute magnitude of this term can be 
given only eifter more data are made available. 

In order to evalmte the term 



the data of Schubauer and Klebanoff were investigated. The magnitude of 
this term is affected by the differentiation and Integration processes 
and for that reason is Influenced by small changes in fairing. A calcu- 
lation was made in which the integrand was believed to be large. The 
calculation indicated that the term 


u2 Jo Jq 



dy dy 


was a small fraction of the term 


1 r 3u*u* 

u2 Jo hx 


cly 


but the resvilts were too inaccurate to warrant presentation. 


In evaluating the term 


U2 



3^u*v' 

3x2 


3y 


the data were sufficient to give a magnitude to but 

3x 


were 


insufficient to permit calculation of 


3x^ 


. Because, however, the 


factor 5 appears in the numerator, it is believed that the magnitude 
of the term is small. 



10 


MCA TN 2571 


CONCUJDIEG REMARKS 


Experimental data show that the streamwise derivative of the turbulent 

longitudinal momentum pu'u’ may he of sufficient magnitude to require 
its inclusion in the application of the Von Karman momentum theorem to 
turbulent boundary layers near separation- The data show that for the 
particular case considered the Von Karman momentum theorem with only 
mean-flow terms gave satisfactory resiilts up to a value of the bou n d a ry - 
layer -velocity -prof lie shape parameter H of approximately 1-5* A 
quantitative discussion of the importance of the fluctuating components 
can be made only when more data become available. 


Langley Aeronautical. Laboratory 

National Advisory Com mi ttee for Aeronautics 
Langley Field, Va., Augiist 29, 1951 
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APPENDIX 

DERIVATION OF THE VON KARMAN MCMENTUM THEOREM 
FOR TURBULENT BOUNDARY LAYERS 


The sin^jlified Navier-Stokea equations, as given previously hy 
equations (4) and (5)^ are 


^ + V ^ 

Sx dyj 


'Su'u' ^u' v' 

4 * — 

^ by 


+ V 

P ^ by^ 


3u' V* ^ ^'v' 


^ ^y 

From the last of these two equations 


1 ^ 
P by 


PO - P8 = P 


Su'v' 


bx 




dy = p 


Su'v' 


<3y 


because for simplicity v'v' is as stoned to approach zero at y = 5 
and at y = 0. Also 


PO 



^'vA 


Then 


P = Ps 


+ p r 

^0 3x 




<3y 


and 


1 Sp 

iP Sx 


3u 3 du'v' 

U — + / 

bx Sx v70 



d [ ^'v' 
5x\ bx 



dy (Al) 


•j 
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"because 

1 ^ = _ u M 

■p ^ Sx 

By su"b 8 tltutlng equation (A 1 ) into the first of the simplified 
Navier-Sto'kes equations, the following equation can he written for a 
turbulent boundary layer in a two-dimensional incompressible flow: 



In order to describe the action of the turbulent boundary layer as a 
whole, equation (A 2 ) is integrated with respect to y through the 
boundary layer: 
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If these quantities are substituted in equation (A 3 ), the following 
relation is obtained: 








ay dy 


[Ah) 


where 


Su'v* 


'0 ^ 


dy = (u’v' )g - (u’v ’ )q = 0 


Integration by parts gives 



Bx 


clyj^ay=u|- 

oy dx 



u dy 


1 a d 5 

2^Jo 


(A5) 
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Substitute equation (A5) in equation (Ah) and divide by U^j then 


-- r 

Sx Jq 


-2. 1 5 



- ^ 1 S6 


1 ^ 

2 u2 Sx 



“^0 _ 1_ 3u’u' ^ S^u'v' 

py2 u 2 Jo ^ Jo Sx^ 


^ I* ay ety 

Jo Jo \ ^x Sy / 

By adding and subtracting 


1 1_ ^ 

2 xj2 ^ 




TT ^ 1 

° - 2 S 


the following equation can be obtained 



(u - U) dy + ^ = 


(a6) 


L. ^ dy - ^ ^ ^ r r ^ 

U2 Jq ax U2 Jq ^^2 Jo Jo ^ / 

(AT) 


Let the momentum loss M of the stream be defined by 


M = pU^e 

where 6 is the momentum thickness, then 
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or 




Now 




(A8) 


where 6*^ the displacement thickness, is defined by 

^ = J (l-^')dy = 5*(x) 


U. 


By using equation (A8), equation (A7) may be written finally as 


“^0 ^ 
pu2 " dx 


(H + 


2 ) ® ® 

U dx 


1 r 8u'u' 

Jo ^ 


dy + 


r ^ fa'v' + 

Xj2 Jo Jo Si\8x 8y / 


dy dy - 


U2 



8x2 


<3y 


where H is defined as 5*/0. This equation is the momentum theorem for 
turbulent boundary layers and is given as equation (6) in the body of 
the paper. 
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Figure 1.- Streamwlse variation of the square or the ratio of the 
velocity at the outer limit of the 'boundary layer to a reference 
velocity. Data vere obtained from reference 11. 
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j Figure 2,- StreamwlBc variation of the boundary- layer momentum thickness. 

j Data vere obtained from reference 11, vo 
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Figure 3»- Streamwlse variation of the houndary- layer shape parameter H. 
Data were obtained from reference 11. 
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Figure 4.- Effect on Ijoundary-layer-Teloclty distribution of the correc- 
tion factor to the pltot-tube measiirementB . Data were obtained from 
reference 11. x = 25.4 feet. 
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Figure 6.- Streamwlse variation of the mean-flow terms in the Von Kajrm^ 
momentum theorem as computed from the data of reference 11. 
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Figure 8. - Streamwlae 
fluctuation ratio 


a nd nor mal variation of the turbulent longitudinal 
Data were obtained from reference 11. 





